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1.  Introduction 


Despite  of  technological  importance  of  the  Si(OOl)  surface  not  much  is  known 
about  an  electronic  structure  of  its  defects.  Hamers  and  Kohler^  studied  atomic-sized 
defects  on  Si(OOl).  Using  the  Scanning  Tunneling  Microscope  (STM),  tunneling 
spectra  i.e.,  current  vs  bias  voltage  dependencies,  were  obtained  for  the  A,  B  and  C 
defects  as  well  as  for  the  regular  Si(001)-(2xl)  surface.  However,  no  attempt  was 
made  to  deconvolute  the  local  density  of  states  (LDOS)  of  the  defects.  Griffith  et 
al?  examined  the  step  geometry  of  Si(00  l)-(2x  1)  by  STM  but  no  electron 
spectroscopy  has  been  accomplished.  An  energy  of  electronic  states  and  a  charge 
distribution  on  Si(OOl)  steps  have  been  calculated  by  Yamaguchi  and  Fujima. 
Unfortunately,  the  difference  between  the  calculated  spectra  and  experimental  data  is 
very  high.'* 

In  this  paper  we  report  on  local  electron  spectroscopy  of  defect  sites  of 
Si(001)-(2xl)  obtained  by  a  new  turmeling  spectroscopy  technique.  Comparative 
Scanning  Tunneling  Spectroscopy,  CSTS.^  This  data  evaluation  technique  reduces 
the  influence  of  the  tip  condition  on  the  sample  LDOS  deconvolution.  We  also 
employ  the  CSTS  to  monitor  variations  of  local  surface  potential  induced  by  defects. 


2.  Experimental 


An  ultrahigh  vacuum  (UHV)  chamber  with  a  background  pressure  of  2x10'^  *  Torr 
was  employed  for  the  experiments.  The  chamber  was  equipped  with  a  STM  (Omicron),  a 
cylindrical  mirror  analyzer  Auger  spectrometer  (Model  15-110,  Physical  Electronics 
Industries  Inc.),  and  a  quadrupole  mass-spectrometer  (lOOC,  UTI). 

Si(OOl)  crystals  (floating  zone  grown,  p-type,  B-doped,  10  and  100  Q-cm,  Virginia 
Semiconductor)  with  dimensions  of  15x4x0.3  mm^  were  cut  from  the  silicon  wafers  by  a 
diamond  scriber,  mounted  on  an  etched-in-HCl  Ta  frame,  ultrasonicated  in  pure  ethanol, 
rinsed  in  deionized  water  and  introduced  into  the  UHV  chamber  through  the  load-lock 
system.  After  overnight  outgassing  at  ~  900  K,  native  oxide  was  removed  from  the 
Si(OOl)  by  a  short  (~  60  sec)  annealing  at  ~  1450  K  followed  by  slow  (~  2  K/sec)  cooling. 
During  annealing  the  background  gas  pressure  was  always  below  lx  10'^  Torr.  STM 
images  revealed  a  well  organized  (2x1)  structure  with  traces  of  Ni  contamination.^ 

The  tips  were  prepared  by  the  electrochemical  etching  of  polyciystalline  tungsten 
wire.  Before  use,  the  tips  were  annealed  in  UHV  by  direct  contact  with  a  tungsten 
filament  preheated  to  ~  2000  K.  The  tip  condition  was  monitored  by  the  STM  imaging 
stability  and  by  the  Si(001)-(2xl)  surface  appearance. 

Negative  sample  bias  potentials  of  1.3-1.4  V  and  tunneling  currents  of  0.25-0.375 
nA  were  used  as  a  reference  for  the  feedback  control  during  spectroscopy.  To  avoid 
silicon  smface  damage  and  to  attenuate  tip  induced  band  bending  ,  the  tunneling  current 
was  always  kept  below  1  nA.  In  order  to  improve  the  signal  to  noise  ratio,  the  data 
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acquisition  time  was  at  least  300  jisec  long  for  the  tunneling  current  measurement  at  each 
tip-sample  bias  voltage.  Every  presented  spectrum  is  averaged  over  more  than  100 
individual  current  vs  voltage  curves  obtained  at  similar  defects. 


The  comparative  LDOS  of  each  type  of  defect  is  calculated  by  the  following  three 
step  procedure  justified  elsewhere^; 

(1)  The  differential  conductivity  vs  bias  voltage  dependence  was  fitted  by  the  asymmetric 
tunneling  probability  function: 

F(S,  V)  =  Aj.t(s,  +  A^t(5  -  (1) 


As  are  the  proportionality  coefficients  related  to  the  tip-surface  effective  contact  area  and 
these  quantities  are  proportional  to  the  tip  and  the  sample  densities  of  states  at  the  Fermi 
level,  respectively;  5  is  the  tip-sample  separation;  m  and  e  are  the  electron  mass  and 

charge,  respectively;  ^  O  average  of  sample  0s  and  tip  0t  work 


functions;  ^  is  the  tuimeling  electron  energy  referenced  to  the  average  of  the  tip  and 
sample  Fermi  potentials;  h  is  Planck’s  constant  divided  by  In.  O ,  At  and  As  are  the 
fitting  parameters.  S  is  assumed  to  be  equal  to  6  A  (see  Section  4). 

(2)  The  differential  conductivity  was  normalized  to  the  asymmetric  tunneling  probability 
function  obtained  fi’om  step  (1).  An  individual  LDOS  is  deconvoluted  at  tiiis  step. 
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(3)  The  difference  between  the  individual  LDOS  determined  for  a  defect  site  and  a 
regular  (2x  1)  site  (the  comparative  LDOS)  was  calculated.  This  step  was  undertaken  in 
order  to  reduce  the  tip  condition  influence  on  the  deconvoluted  LDOS  difference. 

This  procedure  provides  highly  reproducible  comparative  spectra  of  various 
defects  and  minimizes  the  influence  of  the  spectroscopy  parameters  and  the  tip  changes. 

3.  Results 

3.1  Electronic  Characterization  of  Defect  Sites 

An  STM  image  of  a  defective  Si(001)-(2xl)  surface  is  shown  in  Figure  1.  The 
image  was  obtained  in  a  constant  tunneling  current  mode  at  a  positive  sample  bias 
potential  of  1.0  V  and  a  reference  cmrent  of  0.5  nA.  The  majority  of  visible  dimer 
vacancies  are  Ni  induced,  the  “split  off  dimer”  defects.^  In  the  left  top  comer  of  the 
image,  clusters  of  these  “split  off  dimei^’  or  Ni  induced  “vacancy  channel”  defects  are 
seen. 

Two  terraces  of  Si(OOl)  are  shown.  The  A,  B,  and  C  defects'  can  be  found  on 
the  terraces.  Densities  of  the  A  and  C  defects  show  a  general  correlation  with  the 
transient  background  pressure  during  the  crystal  aimealing  and  depend  on  the  UHV 
chamber  condition.  We  believe  that  these  defects  are  induced  by  background  gas 
adsorption  during  the  crystal  cooling  following  its  aimealing. 
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The  B  defect  most  frequently  observed  is  a  part  of  the  “split  off  dimer”  defect. 
Their  density  is  correlated  with  the  Ni  surface  concentration.^ 

The  A  and  B  step  regions*  can  be  seen  along  the  edge  of  the  top  terrace 
(Fig.  1).  The  B  step  is  always  seen  with  the  characteristic  bright  protrusions  on  the 
very  edge. 

Figure  2(a)  presents  the  tunneling  spectra,  cmrent  vs  bias  voltage 
dependencies,  obtained  for  the  10  Q-cm,  p-type,  B-doped,  Si(OOl)  crystal.  Curves 
for  three  different  sirrface  sites  are  shown;  the  regular  (2x1)  dimer,  the  Ni  induced  or 
the  “split  off  dimer”  defect,®  and  the  C  defect.^  Corresponding  deconvoluted 
individual  LDOS  plots  are  presented  in  Figure  2(b).  These  data  for  the  non-defective 
(2x1)  dimer  sites  are  in  reasonable  agreement  with  the  known  LDOS  of  the  ordered 
Si(001)-(2x  1)  surface.'*’^  Comparative  LDOS  for  the  Ni  induced  defect  and  the  C 
defect  are  shown  in  Figure  2(c).  The  spectrum  of  the  regular  (2x  1)  dimer  was  used 
as  a  reference  LDOS. 

In  Figure  3(a)  tunneling  current  vs  sample  bias  voltage  dependencies  for  the 
100  Q-cm,  p-type,  B-doped,  Si(OOl)  crystal  are  presented.  Curves  for  three  different 
surface  sites  are  shown:  the  regular  (2x  1)  dimer,  and  the  A  step  and  the  B  step. 
Corresponding  deconvoluted  individual  LDOS  are  presented  in  Figure  3(b). 
Comparative  LDOS  for  the  A  and  B  steps,  the  Ni  induced  defect  and  the  C  defect  are 
shown  in  Figure  3(c).  In  each  case,  the  spectrum  of  the  regular  (2x1)  dimer  was  used 
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as  a  reference  LDOS.  The  spectrum  of  the  Ni  induced  “vacancy  channel”  defect  is 
similar  to  the  “split  off  dimer”  LDOS  and,  hence,  is  not  shown.  The  A  defect 
spectrum  shows  minor  deviations  from  the  regular  (2x  1)  dimer  LDOS  and  is  also  not 
shown. 

The  best  fit  parameters  obtained  during  differential  conductivity  fitting  by  the 
tunneling  probability  function,  eqn.  (1),  are  presented  in  Tables  1  and  2  for  the  "high"  and 
the  "low"  doped  Si(OOl),  respectively.  If  all  four  fitting  parameters  are  used  (Aj,  As,  S 
and  O)  then  tip-sample  separation  S  and  work  function  O  are  unreasonably  high,  S  ~ 
(13-25)  A  and  <I>~  (9-15)  eV,  respectively.  For  the  feedback  parameters  employed  one 
would  expect  an  5  value  of  ~  (5-6)  A  and  its  minor  variation  (less  than  1  A)  from  one 
surface  site  to  another.  Therefore,  in  order  to  diminish  the  observed  discrepancy  and  to 
constrain  the  fitting  procedure  we  assume  S  =  6  A. 


3.2  Spatial  Characterization  of  C  Defect 

The  cross-section  of  the  C  defect  in  two  orthogonal  directions,  parallel  and 
perpendicular  to  the  dimer  rows  of  the  (2x1)  structure,  are  presented  in  Figures  4(a) 
and  4(b),  respectively.  The  STM  image  was  obtained  with  a  positive  bias  potential  of 
1.0  V  applied  to  the  sample  and  the  reference  feedback  current  of  0.5  nA. 

Depression  with  a  characteristic  lateral  range  from  tiie  C  defect  of  15-20  A  is  seen  in 
the  direction  parallel  to  the  dimer  row.  Fig.  4(a). 
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4,  Discussion 


4.1  Tunneling  Spectra  Fitting  Parameters  for  10  T2-cm  and  100  F2-cm  Si(OOl) 

Tables  1  and  2  present  parameters  which  correspond  to  the  best  fit  of  the 
differential  conductivity  vs  bias  voltage  dependence  by  the  asymmetric  tunneling 
probability  function,  eqn.  (1).  Band  bending  induced  by  the  STM  tip  should  affect  the 
tunneling  current  vs  bias  voltage  dependencies,  especially,  in  the  case  of  low  doped 
silicon.  ’  This  phenomenon  is  not  included  in  equation  (1).  Consequently,  the  fitting 
parameters  (Tables  1  and  2)  and,  to  a  smaller  degree,  the  deconvoluted  LDOS  (Figs.  2 
and  3)  should  be  affected  by  this  factor. 

The  local  surface  potential  or  "local  work  function"  obtained  for  the  "high" 
doped  Si(OOl)  (Table  1)  is  comparable  with  values  reported  by  Hamers  and  Kohler 
(cf.  2.5-3. 0  eV,  ref.  1).  Conversely,  the  "local  work  functions"  obtained  for  the 
"low"  doped  Si(OOl)  crystal  (Table  2)  are  too  low.  Most  likely,  this  discrepancy  is 
due  to  the  strong  tip  induced  band  bending  effect  observed  for  the  100  Q-cm  Si(OOl) 
crystal,  cf.  the  energy  gap  size  in  Figs.  2(b)  and  3(b).  Therefore,  the  tunneling 
probability  fitting  function,  eqn.  (1),  should  be  improved.  The  band  bending  effect 
has  to  be  considered  for  a  low  doped  semiconductor.  Therefore,  in  the  following 
discussion  we  will  focus  on  changes  in  the  local  surface  potential  obtained  for  the 
"high"  doped,  10  Q-cm,  Si(OOl)  ciystal  (Table  1). 

4.2  The  LDOS  Variations  and  tc States 
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Concerning  the  comparative  spectra  of  the  Si(OOl)  defect  sites  (Figs.  2  and  3), 
one  may  emphasize  that  changes  in  LDOS  occur  at  energies  close  to  tc  electronic 
states  of  the  dimerized  surface  Si  atoms.  These  states  are:  the  unoccupied  tC  state  at 
~  (0. 5-1.0)  eV  above  the  Fermi  level,  and  the  occupied  %  state  at  ~  (0. 5-1.5)  eV 
below  the  Fermi  level.'*’^^  Therefore,  the  density  of  these  states  (closest  to  the  Fermi 
level)  is  sensitive  to  the  bonding  geometry  of  surface  Si  atoms  and  to  their  reactions 
with  adsorbates. 

4.3  Ni  Induced  and  "Vacancy  Channel"  Defects 

The  LDOS  of  the  Ni  induced  “split  off  dimer”  and  the  “vacancy  charmel” 
defects  (not  shown)  are  similar.  This  is  in  agreement  with  their  common  ongm.  ’ 
Despite  the  significant  rearrangements  of  surface  Si  atoms  induced  by  the  underlying 
Ni  impurity,^’^^  the  LDOS  experiences  almost  no  changes  with  respect  to  the  normal 
(2x1)  surface.  Nevertheless,  the  local  surface  potential  seems  to  be  reduced  in  the 
vicinity  of  the  “split  off  dimer”  and  “vacancy  channel”  defects  (Table  1).  This 
reduction  in  the  local  surface  potential  is  also  seen  in  the  STM  images  as  a  slight 
elevation  of  the  (2x  1)  surface  around  the  Ni  induced  defects  and,  especially,  in  the 
"vacancy  channel"  vicinity  (Fig.  1).  As  was  shown  by  Koke  and  Monch,*'*  breaking 
of  an  asymmetric  dimer  with  a  large  dipole  moment  may  account  for  a  significant 
local  potential  reduction  of  ~  0.79  eV. 
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4.4  The  C Defect  and  The  B  Steps  ofSi(OOl) 


Spectra  of  the  C  defect,  Figs.  2(c)  and  3(c),  show  a  pronounced  peak  at 
~  0.75  eV  above  the  Fermi  level  and  a  modulation  of  the  density  of  occupied  states 
around  -  1.0  eV  below  the  Fermi  level.  Despite  much  discussion  in  the 
Uterature,*’^®’^^’*^  the  origin  of  the  C  defect  on  Si(001)-(2xl)  is  still  unclear.  It  was 
shown  that  atomic  hydrogen, water*  ^  and  chlorine*^  adsorption  may  cause  the 
appearance  of  the  C  defect.  Unfortunately,  we  are  not  able  to  clarify  the  question 
since  all  above  mentioned  gases  are  present  in  the  background  gas  of  our  UHV 
chamber  during  the  crystal  annealing  and  cooling. 

In  agreement  with  Hamers  and  Kohler,*  we  found  a  reduction  in  the  energy 
gap  size  of  the  C  defect  spectrum.  Fig.  2(a).  Surprisingly  neither  the  height  elevation 
in  the  STM  images  (Fig.  1)  or  a  decay  in  the  local  surface  potential  (Table  1)  were 
observed  around  the  C  defect  (cf.  ref.  1).  Just  the  opposite  was  observed;  the  (2x1) 
dimers  are  slightly  depressed  in  the  C  defect  vicinity  as  seen  in  a  cross-section 
parallel  to  the  dimer  row.  Fig.  4(a).  No  depression  was  observed  in  the  cross-section 
perpendicular  to  the  dimer  row.  Fig.  4(b);  also  no  change  in  local  surface  potential 
was  calculated  (Table  1).  The  difference  in  the  (2x1)  surface  appearance  in  two 
orthogonal  directions  argues  against  the  local  potential  change  or  the  local  defect 
charging  but  rather  suggests  that  the  observed  depression  (or  elevation  ref.  1)  in  the 
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STM  images  due  to  the  variations  in  LDOS  at  specific  biases  in  the  vicinity  of  the  C 
defect. 

Another  possible  reason  for  the  discrepancy  with  the  data  of  Hamers  and 
Kohler  may  be  the  different  dopiog  of  the  Si(OOl)  crystals  used  in  these  two  studies. 

An  apparent  difference  between  the  C  defect  LDOS  obtained  in  this  paper  and 
the  one  observed  by  Boland^®  for  the  H  induced  “bright  ball-like”  sites,  probably, 
rules  out  the  hydrogen  adsorption  as  an  origin  of  the  C  defect  in  our  case.  However, 
it  should  be  mentioned  here  that  an  attempt  to  reproduce  Boland’s  tunneling  spectra 
acquisition  conditions  (feedback  bias  potential  of +0.8  V,  feedback  current  of 
0.15  nA,  bias  ramp  range  of  +  2.5  V)  led  us  to  the  very  high,  ~  30  nA,  tunneling 
current  at  the  highest  and  the  lowest  biases.  At  such  high  tuimeling  current,  damage 
of  the  Si(OOl)  surface  was  regularly  observed  and  the  corresponding  spectra  were 
severely  affected  by  silicon  atom  movement. 

Interesting  features  were  observed  for  the  B  step  LDOS  (Fig.  3).  Overall  the  B 
step  and  the  C  defect  spectra  are  quite  similar.  The  B  step  LDOS  reveals  a  reduction 
in  the  energy  gap  size  and  two  peaks  at  ~  0.75  eV  above  and  below  the  Fermi  level. 
Fig.  3(c).  As  was  mentioned  before,  the  characteristic  bright  protrusions  which  are 
visible  on  the  very  edge  of  the  top  terrace  (Fig.  1)  were  always  observed  during  die 
STM  imaging  of  Si(001)-(2xl),  and  are  umelated  to  the  density  of  the  C  defects  on 
terraces.  Two  explanations  for  the  similarity  of  the  B  step  and  the  C  defect  LDOS 
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may  be  advanced:  (1)  The  C  defects  are  mobile  during  sample  cooling  and 
preferentially  localize  at  the  B  step;  (2)  The  B  step  and  the  C  defect  have  similar 
bond  geometry  and  electronic  structure.  We  favor  explanation  (1)  which  suggests  the 
role  of  an  impurity  adsorbate. 

The  selective  location  of  the  defects  on  the  B  step  edge  is  a  sign  of  special 
chemical  properties  of  these  sites  and  may  explain  preferential  evaporation  and 
chemical  etching  of  the  B  step. 

5.  Conclusion 

A  study  of  local  electron  spectra  of  the  defect  sites  on  the  Si(001)-(2xl)  surface 
was  undertaken  by  employing  Comparative  Scanning  Tunneling  Spectroscopy  (CSTS). 

Ni  induced  defects,  the  A,  B,  C  defects,  and  the  A  and  B  steps  were  studied  (Figs.  2  and 

3). 

The  tip  induced  band  bending  significantly  affects  tunneling  spectra,  especially,  in 
the  case  of  low  doped  silicon.  This  is  evident  firom  LDOS  studied  (Figs.  2  and  3)  and 
from  comparison  of  the  local  surface  potential  (Tables  1  and  2)  for  the  “high”  and  “low” 
doped  Si(OOl)  crystals.  Therefore,  the  tunneling  probability  fitting  function,  eqn.  (1), 
should  be  improved.  The  band  bending  effect  has  to  be  considered  in  the  case  of  low 
doped  semiconductors. 
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A  reduction  of  local  surface  potential  was  found  around  the  Ni  induced  “split  off 
dimer”  and  the  "vacancy  channel"  defects  (Fig.  1  and  Table  1).  The  breaking  of  the 
asymmetric  dimer,  which  has  a  significant  dipole  moment,  may  cause  this  local  potential 
reduction. 

In  agreement  with  Hamers  and  Kohler,*  we  found  a  reduction  in  the  surface  states 
energy  gap  at  the  C  defect,  Fig.  2(a).  At  the  same  time,  neither  elevation  in  the  STM 
images  (Fig.  1)  nor  a  decay  in  the  local  surface  potential  (Table  1)  were  observed  around 
the  C  defect  (cf.  ref.  1).  The  observed  anisotropy  of  the  (2x1)  surface  depression  around 
the  C  defect.  Figs.  4(a)  and  4(b),  suggests  that  the  effect  is  caused  not  by  the  local 
potential  change  or  the  defect  charging  but  rather  by  LDOS  variation  around  the  C  defect. 

The  characteristic  bright  protrusions  observed  on  the  B  step  edge  (Fig.  1) 
reveal  a  spectrum  similar  to  the  C  defect  LDOS,  Fig.  3(c).  Taking  into  account  this 
similarity  we  assume  a  common  origin  of  these  two  defects.  The  selective  location  of 
the  defects  on  the  B  step  edge  is  a  sign  of  special  chemical  properties  of  these  steps 
and  may  explain  their  preferential  evaporation  and  chemical  etching. 
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Figure  Captions 


Figure  1.  400x400  A  STM  image  of  Si(OOl).  The  image  was  obtained  with  a  positive 
bias  potential  of  1.0  V  apphed  to  the  sample.  The  raw  data  without  any  correction  are 
presented.  The  various  defects  are  indicated  by  arrows. 

Figure  2.  Site  specific  electron  spectroscopy  of  Si(OOl)  crystal  (p-type,  B-doped,  10 
Q-cm).  The  spectra  were  obtained  with  a  negative  reference  bias  potential  of  -1.3  V 
applied  to  the  sample  and  a  reference  tunneling  current  of  0.25  nA.  (a)  Tunneling  current 
vs  sample  bias  potential  dependencies,  (b)  Deconvoluted  LDOS.  (c)  Deconvoluted 
comparative  spectra.  The  LDOS  of  the  normal  (2x1)  dimer  was  used  as  a  reference. 

Figures.  Site  specific  electron  spectroscopy  of  Si(OOl)  crystal  (p-type,  B-doped,  100 
Q'cm).  The  spectra  were  obtained  with  a  negative  reference  bias  potential  of  -1.4  V 
applied  to  the  sample  and  a  reference  tunneling  current  of  0.375  nA.  (a)  Tunneling 
current  vs  sample  bias  potential  dependencies,  (b)  Deconvoluted  LDOS.  (c) 
Deconvoluted  comparative  spectra.  The  LDOS  of  the  normal  (2x1)  dimer  was  used  as  a 
reference. 

Figure  4.  Cross-section  of  the  C  defect  on  Si(001)-(2xl):  (a)  parallel  to  the  dimer  row  of 
the  (2x1)  structure;  (b)  perpendicular  to  the  dimer  row  of  the  (2x1)  structure.  Sample 
bias  potential  +1.0  V,  reference  current  0.5  nA. 
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Table  1.  The  best  fit  parameters  obtained  during  the  differential  conductivity  fitting  by 
the  tunneling  probability  function,  eqn.  (1).  The  “high”  (10  Q*cm)  doped  Si(OOl).  The  tip 
-  sample  distance,  S,  was  fixed  to  6A. 


Parameters 

(2x1)  Dimer 

“Split  off  Dimer” 

C  Defect 

Ap 

531 

133 

2819 

As 

1173 

86 

2237 

s(A) 

6 

6 

6 

O  (eV) 

2.3 

1.4 

2.5 
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Table  2.  The  best  fit  parameters  obtained  during  the  differential  conductivity  fitting  by 
the  tunneling  probability  function,  eqn.  (1).  The  “low”  (100  Q-cm)  doped  Si(OOl).  The  tip 
-  sample  distance,  S,  was  fixed  to  6A. 


Parameters 

(2x1)  Dimer 

“Split  off  Dimer” 

C  Defect 

A  Step 

B  Step 

At 

6.8 

21.4 

25.3 

10.6 

17.6 

As 

29.3 

35.1 

47.9 

36.3 

33.5 

s(A) 

6 

6 

6 

6 

6 

O  (eV) 

0.95 

0.99 

1.01 

0.97 

1.02 
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